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EVALUATION  OF  1 5TH -ORDER  HARMONICS  IN  THE  GEOPOTENTIAL 
FROM  ANALYSIS  OF  RESONANT  ORBITS 

by 

D.  G.  King-Hele 
Doreeti  M.  C.  Walker 


SUMMARY 

— Satellite  orbits  contracting  under  the  influence  of  air  drag  experience 
15th-order  resonance  when  the  track  over  the  Earth  repeats  after  15  revolutions. 
If  the  orbital  decay  rate  is  slow  enough,  an  orbit  passing  through  the  resonance 
is  appreciably  perturbed  by  the  effects  of  15th-order  harmonics  in  the  geo¬ 
potential.  We  have  used  the  observed  perturbations  in  23  resonant  orbits,  at 
various  inclinations  to  the  equator,  to  determine  the  harmonic  coefficients  of 
order  15  and  degree  15,  16,  17,  ...  35.  Analysis  of  the  changes  inlorbital 
inclination  on  the  23  orbits  gives  the  harmonics  of  odd  degree,  while  those  of 
even  degree  are  found  from  the  changes  in| eccentr ici ty  on  16  of  the  orbits. 

The  values  derived  are  given  in  Tables  6  and  8.  The  coefficients  of  degrees  15, 
16,  17,  ...  23,  should  be  more  accurate  than  any  previously  obtained;  their 
average  sd  isij.4  y  10“^,  equivalent  to  1  cm  in  geoid  height. 

I '  H-/'' x  i  ..  ‘  * 

Comparisons  with  comprehensive  Earth  models  showed  the  Goddard  Earth  Model 
I0B  to  be  the  best,  and  a  standard  deviation  of  about  3  *  10-t^  in  the  GEM  I0B 
l5th-order  coefficients  is  indicated.  ^  -  -v 
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1  INTRODUCTION 

The  gravitational  potential  of  the  Earth  is  usually  expressed  as  a  double  infinite 
series  of  tesseral  harmonics  depending  on  latitude  and  longitude.  The  order  m  of  the 
harmonics  expresses  the  variation  with  longitude,  and  a  harmonic  of  order  m  has 
m  sinusoidal  oscillations  over  360°  of  longitude.  The  degree  2  of  the  harmonic 
(where  l  >  m)  governs  variations  with  latitude;  these  are  more  complex  and  do  not 
concern  us  here. 

If  the  orbital  period  of  a  satellite  is  such  that  its  successive  ground  tracks  over 
the  Earth  are  360°/m  apart,  so  that  the  track  repeats  after  m  revolutions,  the 
satellite  exhibits  mth-order  resonance  and  the  perturbations  due  to  harmonics  of  order  m 
will  build  up  day  after  day  to  produce  quite  a  large  change  in  some  of  the  orbital 
elements.  This  change  can  be  analysed  to  determine  a  lumped  harmonic  of  order  m  ,  that 
is  a  linear  sum  of  individual  harmonics  of  order  m  and  degree  2^,  2^+2,  2^  +  4,  ..., 
where  2^  =  m  or  m  +  I  (depending  on  the  orbital  element  being  analysed,  and  whether 
m  is  odd  or  even).  By  obtaining  values  of  lumped  harmonics  for  many  resonant  satellites 
at  different  inclinations  to  the  equator,  it  is  possible  to  solve  for  the  individual 
harmonics.  That  is  the  aim  of  this  paper  for  harmonics  of  order  m  =  15  ,  and  the  values 
determined  here  supersede  those  obtained  in  Refs  1  and  2. 

A  satellite  experiencing  15th-order  resonance  has  an  average  height  between  470  km 
(for  near-equatorial  orbits)  and  600  km  (for  inclination  120°),  and  at  these  heights  the 
effects  of  atmospheric  drag  are  appreciable.  So  the  contraction  of  the  orbit  under  the 
influence  of  air  drag  brings  it  to  resonance  and  slowly  draws  it  through  resonance.  The 
lower  the  drag,  the  longer  the  resonance  acts,  and  the  better  the  orbit  is  for  analysis. 
We  have  analysed  23  orbits:  the  longest  resonance  lasts  for  5  years;  but  at  some  inclina¬ 
tions  there  are  no  good  specimens  and  we  have  to  utilise  resonances  that  are  effective 
for  only  about  2  months. 

The  theory  of  the  15th-order  resonance  is  given  in  section  2.  The  analyses  of  the 
23  orbits  are  described  in  section  3.  The  results  for  odd-degree  harmonics  of  order  15, 
from  analysis  of  inclination,  are  presented  in  section  4.3;  and  results  for  even-degree 
harmonics,  from  analysis  of  eccentricity,  or  inclination  and  eccentricity  combined, 
for  16  of  the  23  satellites,  are  recorded  in  section  4.4. 

2  THEORY 

2. I  General  B:a  resonance 

The  longitude -dependent  part  of  the  geopotential  at  an  exterior  point  (r,0,X)  can 

3 

be  written  in  normalized  form  as 


oo  1/ 

f  I  Yffl 

2=2  m=l 


_m  , 

(cos 


[C„  cos  raX  +  S.  sin  mX>N. 
(2m  2m  f  2m 


where  r  is  the  distance  from  the  Earth's  centre,  9  is  co-latitude,  X  is  longitude 

3  2 

(positive  to  the  east),  u  ' a  the  gravitational  constant  for  the  Earth  (398600  km  /a  ) 
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and  R  is  the  Earth's  equatorial  radius  (6378.1  km).  The  P™(cos  6)  are  the  associated 

Legendre  functions  of  order  m  and  degree  l  ,  and  C.  and  S„  are  the  normalized 

im  Urn 

tesseral  harmonic  coefficients:  only  those  of  order  m  *  15  concern  us  here.  The 

3 

normalizing  factor  is  given  by 

„2  -  2(2*  4  DU  -  m)!  ... 

fm  It  +  m)!  ‘  U) 


The  rate  of  change  of  inclination  i  caused  by  a  relevant  pair  of  geopotential 

A, 5 

ItC 


coefficients,  C5m  and  ,  near  B:a  resonance  may  be  written 


3r  ■  ilihr  (f)  W>q(k  "•  1  '  -  jstm)  -  q»>|J  ,  (3) 

where  F.  is  Allan's  normalized  inclination  function^,  G„  is  a  function  of 
imp  >  ipq 

eccentricity  e  for  which  explicit  forms  have  been  derived  by  Gooding'’,  Si  denotes 
'real  part  of'  and  j  »  /-T  .  The  resonance  angle  G>  is  defined  by  the  equation 


*  =  a(oj  +  M)  +  B(ft  -  v)  , 


(A) 


where  to  is  the  argument  of  petigee,  M  the  mean  anomaly,  ft  the  right  ascension  of 
the  node  and  v  the  sidereal  angle.  The  indices  y,  q>  k  and  p  in  equation  (3)  are 
integers ,  with  y  taking  the  values  1,  2,  3,  ...  and  q  the  values  0,  ±1,  ±2,  ...; 
the  equations  linking  i,  m,  k  and  p  are:  m  =  yB  ;  k  =  ya  -  q  ;  2p  =  i-  k. 

At  B:a  resonance  the  m-suffix  of  a  relevant  (C^  ,S^  )  pair  is  given  uniquely  by 
the  choice  of  y  .  The  values  of  l  to  be  taken  must  be  such  that  l  Z  m  and  (£  -  k) 
is  even.  The  successive  coefficients  which  arise  (for  given  y  and  q)  may  usefully  be 
gathered  together  in  a  lumped  form  and  written  as'’ 

k  v  q»k_ 

"  Sfm  •  (5 

l  l 


q,k 


I 


q,k_ 

Cfm 


•  ,  •«.».  n 

where  l  increases  in  steps  of  2  from  its  minimum  permissible  value  l  ,  and  the  Q^’ 

are  functions  of  inclination  that  can  be  taken  as  constant  for  a  particular  satellite; 

and  *  1  when  l  =  l 

i  o 

The  rate  of  change  of  eccentricity  e  caused  by  the  (£,m)  harmonic  near  8:a 
resonance  can  be  written"’ 


de 

dt 


n(  1  -e2)"^ 


(iff.  G. 

\a /  imp  tpq  e  J 


ji-Tn+1(C 


~  jSj,m)  explj(Y«> 


,  (6) 


with  the  same  definitions  as  for  equation  (3) 
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As  the 


Apq 


functions  are  of  order 


(M 


it  is  usually  found  that,  for 


orbits  with  eccentricity  less  than  0.1,  the  terms  with  (y,q)  «  (1,0)  produce  the  most 
important  resonance  effects  on  the  inclination,  though  the  terms  with  (y,q)  ■»  ( 1 ,± 1 )  also 
have  to  be  taken  into  account  if  e  is  greater  than  about  0.03.  In  the  equation  for 
the  eccentricity,  the  relative  importance  of  the  terms  is  largely  decided  by  the  value  of 


1  2 

—  G4pq|q  ~  i (k  +  q)e  J  ,  which  is  of  order  jek  for  q  =  0,  of  order  jf  for  q  »  ±1  , 
2 

and  of  order  if  e  for  q  -  ±2  .  So,  for  the  eccentricity,  the  strongest  effects  are 
usually  caused  by  the  terms  with  (y,q)  =  (1,1)  and  (1,-1). 


Terms  of  30th  order,  with  y  =  2,  may  also  sometimes  need  to  be  taken  into  account, 
but  their  numerical  values  do  not  concern  us  here  and  will  be  the  subject  of  a  separate 
paper. 


2.2  The  15th-order  resonance 

For  15th-order  resonance  (8  =  15,  a  =  1),  equation  (4)  for  the  resonance  angle  $ 
becomes 

$  «  oj  +  M  +  15(0  -  v)  ,  (7) 


and  at  exact  resonance  $  =  0  . 

The  theoretical  equation  (3)  for  variation  of  inclination  may  be  written  as 


di  _  n  /R\ 
dt  sin  i  \a) 


15 


.0,1 


_0,1 


(15  -  cos  i)F]5  ]5  7jc15  sin  ®  “  S]5  cos  $ 


1 7e  /R\-  i-' »0 

+  —  (15)(l)F16,15,8  S15  sin(0-m)+C]5  cos(*  -  u>) 


+  (15  -  2  cos  i) (gF.,  , 


-1.2 


"1,2 


a/ 1 1 6 , !  5 , 7 ) S 1 5  sin(*  +  u)  +  C15  C08(*  +  u,)' 


+  terms  in 


(jfe) cos .  .  . 

- sin  ^  '  qU)) 


IN7 


(8) 


where  only  the  three  main  terms,  with  (y,q)  =  (1,0),  (1,1)  and  (1,-1),  are  given 
explicitly. 

_q,k  _q,k 

The  three  pairs  of  lumped  coefficients  and  appearing  in  equation  (8) 

may  be  written  in  terms  of  the  individual  geopotential  coefficients  (C^.S^)  as 

•  ,  ,  Q  k  - 

indicated  in  equation  (5).  Explicitly,  with  the  Q|V  expressed  in  terms  of  the  F 

_q,k 

functions,  the  C  are 
m 


.0,1 

C!5 


'15,15 


IUjIIjI  (*\2c 


15,15,7 


a/  17,15 


.Li 


15,15,7 


15 


(9) 


k  k 
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J.° 

'I  5 


C  l9F18,15,9  K\2- 

16,15  =  \a)  CI8, 

16,15,8 


2  IF 


15 


20, 15 


I7F 


16,15,8 


15 


(10) 


f1’2  -  C  -  15F'8, 15,8  /R\2- 

15  16,15  J3|  [a)  C18, 15 


16,15,7 


♦  (*fc 

13F.6,.5,7 


a/  20,15 


(11) 


and  similarly  for  S  ,  on  replacing  C  by  S  throughout. 

For  the  15:1  resonance,  the  theoretical  variation  of  eccentricity  given  by 

_q,k  _q,k 

equation  (6)  may  be  written  in  terms  of  the  same  C  and  S  as 

m  id 


L 


de 

dt 


id) 


15 


eF 


15 


/  0,1  _0,1  \ 
,15,7^15  sin  *  ~  S15  C0S  y 


.1,0 

sin(4  -  a)  +  C]3  cos($  -  m) > 
_-l,2 

sin($  +  u)  +  C]3  cos(i!>  +  w) 


+  terms  xn 


iioj 


|q|- 


(|q|)i 


jq  -  J(k  +  q)e  } 


(12) 


Three  terms  are  given  explicitly  in  equation  (12),  those  with  (y,q)  =  (1,0),  (1,1)  and 
(1,-1).  The  main  terms  are  expected  to  be  those  with  (y,q)  =  (1,1)  and  (1,-1),  but  the 
term  with  (y,q)  »  (1,0)  is  also  given,  for  consistency  with  equation  (8). 

3  THE  ANALYSIS  OF  THE  23  RESONANT  ORBITS 

The  methods  of  analysis  have  been  explained  in  several  previous  papers,  most 
recently  Refs  6  and  7,  and  the  explanations  will  not  be  repeated  here.  Basically,  the 
observational  values  of  inclination  are  cleared  of  irrelevant  perturbations  and  fitted 
using  the  computer  program  THROE  with  an  integrated  form  of  the  theoretical  equation  (8), 
with  extra  terms  when  appropriate,  to  determine  values  of  the  lumped  coefficients. 
Similarly  the  observational  values  of  eccentricity,  cleared  of  perturbations,  are  fitted 
with  an  integrated  form  of  equation  (12),  with  extra  terms  as  necessary.  With  a  few 
satellites  it  is  useful  to  make  a  simultaneous  fitting  of  inclination  and  eccentricity 
using  the  SIMRES  program.  In  making  the  fittings,  we  regard  20  values  of  inclination 
(or  eccentricity)  as  the  minimum  permissible,  and  we  try  to  analyse  orbits  over  a  period 
of  time  when  ♦  lies  between  -10  and  +10  deg/day,  though  on  some  high-drag  orbits 
larger  values  of  $  have  to  be  allowed. 

In  seeking  resonant  satellites  for  analysis,  the  aim  has  been  to  cover  the  widest 
possible  range  of  inclination  and  to  leave  the  smallest  possible  gaps  in  the  coverage. 
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o on  i  o 
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No  suitable  orbits  at  inclinations  less  than  30°  were  found,  but  this  is  not  so  bad  as  it 
seems,  because  such  orbits  are  influenced  primarily  by  harmonics  of  very  high  degree 
(35-55,  or  even  higher  for  orbits  very  near  the  equator)  and  our  evaluations  only  extend 
to  degree  35.  The  analyses  of  the  23  orbits  are  described  in  sections  3.1  to  3.23  in 
order  of  increasing  inclination.  The  values  of  lumped  harmonics  given  in  sections  3.1  to 
3.23  are  used  (in  section  4)  to  evaluate  individual  harmonic  coefficients.  The  values  of 
the  lumped  harmonics  are  used  with  their  standard  deviations  unchanged,  unless  otherwise 
specified.  The  values  used  in  the  solutions  are  listed  in  Tables  1  and  2. 

2 

In  evaluating  the  fittings  we  often  refer  to  the  measure  of  fit,  e  ,  where  e  is 
defined  as  the  sum  of  squares  of  weighted  residuals  divided  by  the  number  of  degrees  of 
freedom.  The  weighted  residual  is  the  residual  of  an  individual  value  (of  i  or  e)  divided 
by  its  assumed  standard  deviation.  For  the  US  Navy  orbits  used  in  many  of  the  fittings, 
the  assumed  standard  deviation  is  0.003°  in  inclination  and  0.00004  in  eccentricity;  for 
the  RAE  orbits  determined  by  PROP,  the  standard  deviation  given  by  PROP  is  used;  for 
other  orbits  the  standard  deviation  is  as  specified  in  the  appropriate  section. 


Pegasus  1,  1 965-09A  (i  =  31.76 


0.007) 


This  satellite  passed  through  exact  15th-order  resonance  on  8  December  1974  and  the 
inclination  and  eccentricity  were  analysed  over  the  period  September  1974  to  May  1975. 
During  this  time  the  rate  of  change  of  the  resonance  angle,  £  ,  increased  from  -5  to 
+5  deg/day. 


There  were  37  US  Navy  orbits  available  over  the  period  to  be  analysed.  The  values 
of  inclination,  cleared  of  all  perturbations  except  those  due  to  resonance,  were  first 
fitted  with  (y,q)  =  (1,0)  only,  the  (l,±l)  terms  not  being  required  as  the  eccentricity 
was  only  0.007.  The  (1,0)  fit  was  very  satisfactory,  e  =  0.351,  and  the  values  of  the 
lumped  harmonics  well  determined.  A  second  fit  was  tried  with  (y,q)  =  (1,0)  and  (2,0), 
but  this  was  not  acceptable  because  the  value  of  e  increased  and  the  C^g  and  S^g 
coefficients  were  indeterminate.  The  (y,q)  =  (1,0)  solution  gave  the  following  values: 


oJM 

10  C)5  =  30980  ±  1960  , 


10  S|5  =  13540  ±  960 


The  values  of  inclination  and  the  fitting  by  THROE  are  shown  in  Fig  1. 

The  values  of  eccentricity  were  fitted  with  (y,q)  =  (1,1)  and  (1,-1),  as  these 
are  the  most  likely  terms  to  be  required  in  fitting  eccentricity  alone.  The  result  was 
very  disappointing  and  none  of  the  terms  was  determined  with  sufficient  accuracy  to 
be  useful,  the  value  of  e  being  4.9.  A  further  run  with  a  (0,1)  term  included  was 
tried.  This  had  been  required  in  previous  analyses^’^  in  order  to  remove  an  oscillation 
from  the  US  Navy  values  of  eccentricity,  which  arises  because  of  a  wrong  interpretation 
of  the  zonal  harmonic  oscillation  removed  by  the  US  Navy  from  their  values.  However, 
this  technique  was  not  successful  here  and  the  values  of  the  lumped  coefficients  were 
not  acceptable. 


Values  of  lumped  harmonics  (C 
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3.2  OSO  6  rocket,  I969-68B  (i  =  32.97°,  e  =  0.004) 

The  new  fitting  of  the  inclination  by  THROE,  with  (y,q)  =  (1,0)  shown  in  Fig  2,  is 

2 

virtually  identical  to  the  previous  fitting  .  The  curve  fits  the  points  very  well,  with 
e  =  0.28,  and  the  values  of  the  lumped  harmonics  are: 

q_°.>  q_°,l 

10  C15  =  20340  t  750  10  S15  =  6280  +  910  • 

A  fitting  with  (y,q)  =  (1,0)  and  (2,0)  was  also  tried,  but  the  (2,0)  terms  were 
indeterminate.  The  value  of  f1  increased  from  -3.7  deg/day  initially  to  17.3  deg/dav 
at  the  end  of  the  analysis. 

In  analysing  the  eccentricity,  a  (y,q)  =  (0,1)  term  was  added  for  the  reason  given 
in  section  3.1.  The  fitting,  with  (y,q)  =  (1,1),  (1,-1)  and  (0,1)  terms,  was  poor,  with 
c  =  3.1.  None  of  the  values  of  lumped  harmonics  was  more  than  twice  its  standard  devia¬ 
tion.  So  the  attempt  was  abandoned. 

3.3  San  Marco  I,  1964-84A  (i  =  37.80°,  e  =  0.042) 

This  is  the  only  15th-order  resonant  orbit  available  at  an  inclination  between  33° 

and  43°,  so  it  had  to  be  utilized  if  possible.  By  an  unfortunate  chance,  the  values  of 
_0,1  _0 , 1 

both  and  are  very  small  at  an  inclination  of  38  ,  so  the  variation  in 

inclination  is  also  very  small.  To  compensate  for  this  ill  luck,  however,  28  accurate 

orbits  at  dates  near  resonance  were  available  from  the  archives  of  the  Smithsonian 

2 

Astrophysical  Observatory,  and  these  were  analysed  in  1974.  The  last  of  the  28  orbits 
does  not  fit  well,  and  better  values  for  the  lumped  harmonics  have  now  been  obtained  by 
omitting  the  28th  orbit.  The  values  of  inclination,  and  the  fitting  by  THROE,  are  shown 
in  Fig  3.  The  value  of  $  runs  from  -27  to  +25  deg/day.  Since  the  eccentricity  is 
appreciable  (0.042),  we  have  to  use  all  three  (v,q)  terms  -  (1,0),  (1,1)  and  (1,-1).  The 
assumed  accuracy  in  the  values  of  i  was  0.001°  and  the  THROE  fitting  gave  £  =  0.50, 
with  the  following  values  for  the  lumped  coefficients: 

q_0  > 1  q_0,l 

10  C|5  =  560  +  580  10  S15  =  "  2000  ±  1450  ' 


The  eccentricity  was  also  analysed,  taking  (y,q)  =  (1,1)  and  (1,-1).  The  assumed 
accuracy  was  0.00002  and  the  value  of  e  was  2.9.  Unfortunately  in  this  fitting  and  in 
a  simultaneous  fitting  of  i  and  e  with  the  SIMRES  program,  the  values  of  lumped 
harmonics  obtained  were  not  accurate  enough  to  be  acceptable. 

3.4  HEAP  3,  1 979-82A  (i  =  43.60°,  e  =  0.001) 

The  third  high-energy  astronomical  observatory  was  launched  on  20  September  1 9 7 Q 
into  an  orbit  very  close  to  l5th-order  resonance.  Exact  resonance  was  reached  on 
14  November  1979.  NASA  orbits  are  available  at  2-day  intervals  and,  although  it  was 
feared  that  orbital  manoeuvres  might  have  disturbed  the  effects  of  the  resonance,  the 
analysis  proved  to  be  quite  satisfactory.  Fig  4  gives  the  values  of  inclination  cleared 
of  perturbations,  the  fitting  by  THROE  with  (y,q)  -  (1,0)  being  shown  as  a  broken  line. 
With  the  accuracy  of  i  taken  as  0.002°,  the  value  of  r  was  0.63.  This  fitting  follows 
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Che  main  trends  of  the  variation,  but  there  is  obviously  an  unmodelled  oscillation  with  a 
period  about  half  that  of  the  argument  of  perigee  oj  .  So  THROE  was  run  with  (y,q)  =  (1,0) 

and  (0,2),  that  is,  with  sin  2ui  and  cos  2io  terms  added.  The  value  of  t.  decreased 

to  0.49,  the  values  of  the  lumped  harmonics  changed  by  less  than  1  sd,  and  their  standard 
deviations  were  reduced  by  about  20%.  This  solution,  shown  by  the  unbroken  line  in  Fig  4, 
was  preferred,  and  gives: 

qJM  qJM 

10  C]5  =  -  467  ±  34  10  S  =  -  767  ±  106  • 

The  values  from  the  (y,q)  =  (1,0)  fitting  were  -504  ±  41  and  -666  ±  132  respectively. 

These  THROE  runs  used  53  values  of  i  ,  with  $  increasing  from  -7.7  to  +6.8  deg/day,  and 

gave  better  results  than  runs  with  55  and  57  values. 

The  values  of  eccentricity  from  the  same  53  orbits,  cleared  of  air-drag  perturba¬ 
tions,  were  also  successfully  fitted  with  (y,q)  =  (1,1)  and  (1,-1),  after  a  residual 
oscillation  correlated  with  w  was  removed  by  using  a  value  of  0.586  for  lO^J^.  The 
assumed  accuracy  was  0.00004  and  e  was  0.39.  Fig  5  shows  the  values  cleared  of 
perturbations  and  the  fitted  curve.  The  values  of  the  lumped  harmonics  are: 


9-'.° 
10C!5  - 

-  860  + 

150 

q_  1  » 0 

10  S,5  = 

-  1930  t  160 

9_-',2 

10  Cl 5  " 

-  234  ± 

34 

9--', 2 

10  Sl5  = 

185+67  . 

Although  the  fitting  of  the  curve  in  Fig  5  is  not  perfect,  only  the  fourth  value  needed 
its  standard  deviation  increased,  by  a  factor  of  4,  when  the  values  were  used  in  the 
solutions  for  individual  coefficients. 

3 . 5  Tournesol  1  rocket,  1971-30B  (i  =  46.36°,  e  =  0,011) 

Toumesol  1  rocket  passed  through  exact  15th-order  resonance  on  5  August  1978  and 
over  a  period  from  May  to  October  there  were  24  US  Navy  orbits  available  for  analysis. 
During  this  time  the  value  of  £  changed  from  -9  to  +10  deg/day. 

The  24  values  of  inclination,  cleared  of  perturbations  except  those  due  to 
resonance  were  fitted  by  THROE  with  (y,q)  =  (1,0).  This  fitting  gave  e  =  0.423  and 
the  values  of  the  lumped  harmonics  were : 

q_0»>  q_0,i 

10  C]5  =  -  596  ±  77  10  SI5  =  "  869  1  47  • 

The  values  of  inclination  and  the  THROE  fitting  with  (1,0)  are  plotted  in  Fig  6,  and 
although  the  fit  looks  good,  both  standard  deviations  had  to  be  increased  by  a  factor  of  4 
in  the  solutions.  Two  further  runs,  with  (y,q)  =  (1,0)  and  (2,0),  and  with  (y,q)  = 

(1,0),  (1,1)  and  (1,-1),  produced  larger  standard  deviations  and  increased  c  . 


81006 


harmonics  (C,S).!  and  (C,S)  ’  for  the  16  satell 


-54 A  90.21  -92  ±  48  -170  t  112*  -0.01237  1.1  '  0.6 

-62.9  £  2.6  -53.4  •  1.6  -0.3833 
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The  values  of  eccentricity  were  fitted  by  THROE  using  (y,q)  =  (1,1)  and  (1,-1),  and 
this  fitting  was  accepted  as  the  addition  of  a  (0,1)  term  gave  no  advantage.  The  values 
of  the  lumped  harmonics  were: 


-  2020  ±  280 


oJ.O 

qj»0 

10C,5  * 

86  i  350 

,0  Sl5 

9_-1 »2 

10C15  = 

-  6  ±  98 

q_-«,2 

.0  s15 

The  standard  deviations  of  C, 


J,  0 

and  S , _ 


175  i  74 


had  to  be  increased  by  factors  of  2  and  4 


respectively.  The  values  of  eccentricity  and  the  THROE  fitting  are  plotted  in  Fig  7. 

A  SIMRES  fitting  was  tried  using  (y,q)  =  (1,0),  (1,1)  and  (1,-1)  for  both  inclina¬ 
tion  and  eccentricity.  This  gave  similar  results  for  eccentricity,  but  doubled  the 
standard  deviations  and  e  in  the  fitting  of  inclination,  so  the  separate  fittings  were 
pre ferred. 

3.6  Intercosmos  11,  I974-34A  (i  =  50.64°,  e  =  0.002) 

This  satellite  passed  through  exact  15th-order  resonance  on  I  October  1976  and  was 
analysed  over  a  two-year  period  by  Walker^.  The  values  for  the  lumped  coefficients  from 
the  fitting  by  THROE  of  inclination  and  eccentricity  were  taken  from  equations  (II)  and 
(16)  of  Ref  7.  They  are  given  in  Tables  1  and  2. 

3.7  Explorer  44  rocket,  1971-58B  (i  =  51.05°,  e  =  0.011) 

.  9 

The  orbit  of  1971-58B  near  15th-order  resonance  was  determined  by  Hiller  from  700 
observations  using  the  RAE  computer  program  PROP.  The  drag  was  rather  high,  so  there 
was  no  chance  of  results  as  good  as  those  from  1974-34A.  Mixed  PROP  and  US  Navy  orbits 
were  used  in  fitting  the  inclination,  and  the  analysis  was  improved  by  subtracting  0.002° 
from  all  the  US  Navy  values  of  inclination  and  increasing  their  standard  deviation  to 
0.005°.  Previously,  we  used  the  fitting  with  (y,q)  =  (1,0),  (2,0),  (1,1)  and  (1,-1),  but 
in  the  light  of  subsequent  experience  we  regard  the  (2,0)  terns  as  dubious,  and  we  have 
chosen  the  solution  with  (y,q)  =  (1,0),  (1,1)  and  (1,-1).  The  fitting  is  shown  in  Fig  8 
and  the  values  of  the  lumped  harmonics  obtained  are: 


354  ±  47 


=  -  248  ±  45 


Since  the  fitting  (c  =  1.13)  is  not  as  good  as  might  be  hoped,  and  since  accurate  values 

are  available  from  1974-34A  at  nearly  the  same  inclination,  some  relaxation  of  the 

_0,1 

standard  deviations  is  reasonable,  and  that  of  was  doubled  in  the  final  solutions. 

Hiller  also  analysed  the  eccentricity  of  this  satellite,  and  his  fitting,  with 
(v,q)  ”  (1,1)  and  (1,-1),  is  shown  in  Fig  9.  Again,  the  determination  of  four 
coefficients  by  fitting  28  points  is  not  likely  to  be  reliable,  but  we  decided  to  use  the 
values  and  to  increase  the  standard  deviation  if  the  residuals  in  the  solutions  were  poor. 
The  original  values  were: 


1006 


Jf#*- 


=  -  466  ±  58 


253  ±  60 


=  -  50  ±  23 


9-~' 

10  S>5 


=  45+14 


In  Che  final  solutions,  the  standard  deviation  was  quadrupled  for  C. 

J,0  .-1,2 

S , r  and  C . ,  . 


and  doubled  for 


3.8  Ariel  1,  1962-15A  (i  =  53,82  ,  e  =  0.022) 

Ariel  1,  the  world's  first  international  satellite,  passed  through  exact  15th-order 
resonance  on  8  May  1973.  The  inclination  and  eccentricity  were  analysed  by  Walker*^*  over 
a  period  of  six  months  centred  on  the  exact  resonance. 

Mixed  PROP  and  US  Navy  orbits  were  used  in  the  fitting  of  the  inclination  and  we 
have  used  the  solution  with  (y,q)  =  (1,0)  and  (1,1).  This  was  the  solution  recommended 
by  Walker  as  best  and  the  values  of  the  lumped  harmonics  are: 


=  -  370  +  14 


-  114  +  31 


The  values  of  lumped  harmonics  obtained  from  the  fitting  of  eccentricity  by  THROE 

were  not  used  in  the  previous  determination  of  the  even  harmonic  coefficients*.  Here  we 

use  the  values  from  a  fitting  with  (y,q)  =  (1,1)  and  (1,-1)  given  in  Run  9  of  Table  5  by 
_ 10  _ _ 


Walker 


The  values  are: 


oJ.O 

10  C,5 


=  -  76  ±  18 


=  172  .+  38 


10'C)5  =  151  ±  15  10  S  =  11+34  . 

(These  values  are  obtained  from  the  coefficients  -B,  A,  -D  and  C  respectively  in  Ref  10 

after  dividing  by  -0.7602  for  the  first  pair  and  1.581  for  the  second  pair.)  In  the 

solutions  for  the  individual  harmonics  it  was  foui  d  necessary  to  increase  the  standard 

_  1,2  _1,0 

deviations  of  and  by  factors  of  4  and  2  respectively. 

3.9  Cosmos  72,  1965-53B  (i  =  56.04°,  e  =  0.003 


This  orbit  passed  slowly  through  resonance  during  1972,  with  a  change  in  inclina¬ 
tion  of  0.07°.  The  previous  analysis  by  THROE,  using  seven  PROP  orbits**  and  45  Navy 
orbits  with  (y,q)  =  (1,0),  gave  the  following  values  of  the  lumped  harmonics: 


10  C  =  -  233.4  +  3.3 


=  -  103.4  t  8.4  . 


The  fit,  which  is  excellent,  with  t  =  0.7,  is  shown  in  Fig  6  of  Ref  2.  These  values 
fitted  well  in  the  previous  solutions  and  the  C  value  was  the  most  accurate  previously 
obtained.  However,  we  found  that  the  S  coefficient  did  not  fit  our  new  solutions, 
which  require  a  value  near  -60  rather  than  -103.  So  further  fittings  with  THROE  were 
tried,  including  extra  terms  such  as  (2,0),  (1,1)  and  (1,-1);  omitting  the  last  6  and 
then  the  last  12  values;  and  then  adding  0.0005°  to  the  US  Navy  values.  The  values  of  the 
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lumped  coefficients  obtained  did  not  differ  significantly  from  those  quoted  above,  and  the 
numerical  value  of  the  S  coefficient  never  fell  below  101.  So  the  old  values  were 
retained,  but  the  standard  deviation  of  the  S  coefficient  had  to  be  increased  by  a 
factor  of  4.  The  non-conforming  S  value  is  puzzling,  and  we  plan  to  determine  PROP 
orbits  throughout  the  resonance  phase  in  the  hope  of  resolving  the  problem. 

In  fitting  the  values  of  eccentricity,  the  mismatch  between  PROP  and  US  Navy  values 
was  removed  by  subtracting  0.0001  sin  u>  from  the  PROP  values.  The  fitting  with  (y,q)  » 
(1,1)  and  (1,-1)  was  fairly  satisfactory,  with  e  =  1.62,  but  as  expected  there  was  an 
unmodelled  variation  with  the  same  period  as  ui  .  The  addition  of  (y,q)  *  (0,1)  terms 
led  to  a  much  improved  fit,  with  e  =  0.80,  shown  in  Fig  10.  The  values  of  the  lumped 
harmonics  were  as  follows: 


9J.O 

9J.° 

10  C15 

=  18+17 

10  S15 

=  57  ±  23 

9-'1’2 
10  C15 

=  106.9  ±  8.7 

io  s,5 

=  2.4  ±  8. 

3.10  Tiros  7  rocket,  1963 

-24 B  (i  =  58.20°, 

e  =  0.002) 

This  satellite,  like  1974-34A,  was  analysed  by  Walker2  over  a  two-year  period,  exact 
15th-ordir  resonance  occurring  on  3  March  1977.  The  lumped  coefficients  from  the  THROE 
fittings  of  inclination  and  eccentricity  are  taken  from  equations  (9)  and  (14)  of  Ref  7. 
The  values  are  given  in  Tables  1  and  2. 

3.11  Cosmos  373,  1970-87A  (i  =  62.92°,  e  =  0.007) 

The  orbit  of  Cosmos  373  was  already  past  resonance  when  its  initial  manoeuvres 
ceased;  apart  from  this  defect,  it  is  a  good  satellite  because  its  decay  rate  was  slow. 
Since  there  are  no  ideal  resonant  orbits  at  inclinations  between  58.2°  and  74.0°,  we  have 
to  use  four  imperfect  specimens,  of  which  this  is  the  first. 

.  .  .  .  2 

In  1974  we  analysed  the  variations  in  inclination  using  24  US  Navy  orbits  .  The 

work  showed  that  only  (y,q)  »  (1,0)  could  be  used;  the  addition  of  (1,1)  and  (1,-1)  terms 
was  disastrous  because  of  the  correlations  caused  by  the  near-constancy  of  w  .  A  new 
fitting  of  the  orbits  with  (y,q)  *  (1,0)  has  been  made,  using  improved  methods  for 
removing  lunisolar  and  air  drag  perturbations.  The  worst-fitting  point  was  the  last, 
so  it  was  omitted.  Fig  1 1  shows  the  curve  fitted  by  THROE  to  the  23  orbits,  with  $ 
increasing  from  0.7  to  4.6  deg/day.  The  fitting  is  very  good,  with  t  =  0.31,  and  gave 
the  following  values  for  the  lumped  coefficients: 

q_0,l  a_0,l 

10  Cj5  =  -  5.3  t  3.2  10  SI5  =  “  32-8  ±  2 • 5  . 

The  eccentricity  was  analysed  in  1974  but  the  fitting  of  the  curve  was  very  poor 
(see  Fig  5  of  Ref  1),  and  we  decided  it  was  not  now  acceptable. 

1 2 

The  orbit  of  1970-87A  has  been  determined  by  Brookes  at  selected  epochs  between 
1970  and  1975,  and  the  first  five  of  his  epochs  fall  within  the  time  interval  of  our 
analysis.  Unfortunately,  it  was  not  possible  to  mix  these  orbits  with  the  US  Navy  orbits. 


J 
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and  further  orbits  are  now  being  determined  at  the  University  of  Aston  in  the  hope  of 
defining  the  variations  in  inclination  and  eccentricity  more  precisely. 

3.12  Tansei  3  rocket,  1977—1 2 B  (i  =  65. 49°,  e  =  0.029) 

When  the  Japanese  Tansei  3  satellite  was  launched  on  19  February  1977,  a  rocket, 
1977-12B,  was  left  in  a  lower  orbit.  The  rocket  passed  through  15th-order  resonance  on 

16  March  1978  and  decayed  on  21  March  1979.  Because  of  its  high  drag,  this  orbit  is  far 

from  ideal  for  resonance  analysis. 

In  the  THROE  fitting,  25  weekly  US  Navy  orbits  were  used,  covering  a  range  of  t 
so  wide  (-19  deg/day  to  +29  deg/day)  that  several  oscillations  in  the  'tail'  of  the 
resonance  are  inevitably  included.  Fig  12  shows  the  curve  fitted  by  THROE.  As 

expected,  the  change  at  resonance  was  small,  but  the  fitting  was  good  (e  =  0.54)  and 

the  values  of  the  lumped  harmonics  were: 

q-O.'  q  0, ! 

10yC15  =  -  34  .+  7  10  si5  =  -  18  r  14  . 

In  the  solutions  for  the  individual  coefficients,  the  standard  deviation  of  the  first  of 
these  values  had  to  be  doubled  but  the  second  fitted  well. 

Analysing  the  variation  in  eccentricity  seemed  sure  to  be  fruitless  and  was  not 
attempted . 

The  orbit  of  1977-1 2B  is  now  being  determined  with  PROP  at  the  University  of  Aston 
and,  when  this  work  is  completed,  more  accurate  results  should  be  obtained. 

3.13  Cosmos  462,  1971-106A  (i  =  65.70°,  e  =  0.045) 

The  inclination  and  eccentricity  of  this  high-drag  satellite  were  analysed  by 
1  3 

Walker  and  the  results  from  that  analysis  are  used  here.  The  inclination  was  fitted  ' 
THROE  with  (y  ,q)  =  (1,0)  and  the  values  obtained  were: 

9_°+ 1  9-°» ' 

10X(5  =  -  36  *  21  10  S15  =  9  i  17  • 

The  eccentricity  values  were  also  analysed  and  the  SIMRES  fitting,  using  (y,q)  - 

1  3 

(1,1)  and  (1,-1),  was  recommended  by  Walker  as  the  best  solution.  The  values  were: 

g_' >0  q_* 

10  C|3  =  51-24  10  SI5  =  -55+7 

-12  - i  2 

109C|5’  =  -  67  +10  1 09S  j  5  =  -  19  .+  24  . 

_  1,2  _1 ,0 

For  the  C.„  and  SIC  values  it  was  necessary  to  increase  the  standard  deviations 
li  15 

by  factors  of  4  and  10  respectively. 

3.14  China  2  rocket,  197I-I8B  (i  =  69.84°,  e  °  0.040) 

China  2  rocket  was  in  orbit  from  3  March  1971  until  16  February  1976.  Its  orbit 

1  4 

has  been  determined  by  Hiller  at  114  epochs  from  moie  than  7000  observations,  using 
the  PROP  orbit  refinement  program.  Hiller  analysed  the  orbital  changes  at  four  resonances. 
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and  his  analysis  of  inclination  at  15th-order  resonance  with  (y,q)  “  (1,0)  gave  the 
following  values  of  lumped  harmonics: 

o_0, 1  g.O.l 

10  C|5  -  -  37  ±  6  10  S15  *  10  i  6  • 

The  fitting  is  shown  in  Fig  11  of  Ref  14.  Since  this  was  a  high-drag  orbit,  some  relaxa¬ 
tion  of  the  standard  deviations  is  likely  to  be  necessary;  in  fact  the  second  was  doubled, 
The  values  were  still  most  useful,  however,  since  this  is  our  only  satellite  at  inclina¬ 
tions  between  65.8°  and  74.0°. 

Hiller  also  attempted  to  analyse  the  variations  in  eccentricity,  but  the  resulting 

14 

values  of  lumped  coefficients  were  indeterminate 

3.15  Cosmos  387,  1 970-1 1 1 A  (i  -  74.00°,  e  =  0,001) 

This  low-drag  satellite  gave  excellent  results from  analysis  of  19  PROP  orbits 
and  55  US  Navy  orbits  between  May  1971  and  July  1972.  Here  we  use  the  same  values  of  the 
lumped  coefficients  as  before: 


Q.0,1 

n_0,l 

.0  c15  = 

-  26.0 

± 

1.0 

10  S15  - 

-5.2±l 

1.3 

qj.0 

-  18.0 

oJ»0 

1°  C.5  - 

± 

3.3 

10  S>5  " 

-  44.1  + 

2.5 

9_->,2 

o_-l,2 

10  C.5  - 

-  46.5 

± 

2.7 

,0S>5  " 

-  40.5  ± 

4.0 

The  first  two  coefficients  come  from  a  fitting  of  inclination  with  (y,q)  =  (1,0)  and 

(2,0).  The  last  four  coefficients  come  from  a  fitting  of  eccentricity,  with  (v,q)  =  (1,1) 

and  (1,-1),  being  obtained  from  the  coefficients  -B,  A,  -D  and  C  respectively  in  Ref  1, 

after  dividing  by  -0.6116  for  the  first  pair  and  by  0.4001  for  the  second  pair.  The 

_1,0 

standard  deviation  of  had  to  be  increased  by  a  factor  of  10;  the  reason  for  this 

discrepancy  is  not  known,  but  may  be  connected  with  the  very  unusual  variation  of  e 
resulting  from  its  very  low  value  (see  Fig  2  of  Ref  15). 

3.16  Cosmos  395  rocket,  1971-13B  (i  =  74,05°,  e  =  0.002) 

This  satellite  is  almost  a  twin  of  1 970— 1 1 1 A  and  the  results,  obtained'^  from 
analysis  of  21  PROP  orbits  and  67  US  Navy  orbits  between  September  1971  and  October  1972, 
were  even  better  than  for  1970-1 1 1A.  We  used  the  same  values  of  the  lumped  coefficients 
as  before: 


q_0 , 1 

,0C>5  = 

-  24.6  +  1.3 

q_0,l 

1°  s15 

=  ■ 

-  6.1  ± 

1  .0 

oJ.O 

,0C15  = 

-  19.8  +  1 .8 

oJ.O 

10  S15 

=  ■ 

-  24.8  ± 

0.7 

9--'. 2 
,0C15  = 

-  45.5  i  2.0 

9-"'. 2 

.o  s15 

=  ■ 

-  35.2  + 

1  .0 

• 

The  first  two  coefficients 

come  from  a  fitting 

of  inclination 

with  (y 

.q)  = 

0,0) 

and 

(2,0).  The  last  four  are 

from  a  fitting  of  eccentricity 

with 

(Y,q)  - 

(1,1) 

and 

(1,-0, 
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and  are  obtained  from  the  coefficients  -B,  A,  -D  and  C  respectively  in  Ref  1,  after 
dividing  by  -0.6104  for  the  first  pair  and  by  0.3976  for  the  second  pair. 

3.17  Cosmos  956  rocket,  1977-95B  (i  =  75.82°,  e  =  0.029) 

This  is  the  most  recent  resonance  to  be  analysed:  exact  15th-order  resonance 
occurred  on  19  May  1980.  Both  the  inclination  and  eccentricity  have  been  analysed  over  a 
period  from  mid  March  until  the  end  of  August  1980  using  24  US  Navy  orbits;  during  this 
time  $  changed  from  -11  to  +12  deg/day. 

The  values  of  inclination  cleared  of  non-resonant  perturbations  are  plotted  in 
Fig  13  and  the  curve  shows  the  THROE  fitting  to  the  values  with  (y,q)  «  (1,0),  (1,1)  and 
(1,-1),  the  (1,1)  and  (1,-1)  terms  being  included  as  the  eccentricity  is  0.029.  The 
lumped  coefficients  are  as  follows: 

o_0J  9_°. 1 

.°C15  -  -22.5+5.1  10S15  -  -3.0  ±5.4  . 

The  values  of  eccentricity  were  fitted  with  THROE  over  the  same  period  with  (y,q)  = 
(1,1)  and  (1,-1).  Adding  a  (0,1)  term  was  tried  for  the  reason  given  in  section  3.1  but 
the  values  were  indeterminate.  Another  fitting  with  the  (1,0)  term  included  was  made 
and  used  with  the  THROE  fitting  of  inclination  in  a  SIMRES  fitting  of  inclination  and 
eccentricity  together.  This  procedure  proved  successful,  as  the  addition  of  the  e-terms 
from  equation  (8)  helped  to  provide  a  better  determination  of  the  lumped  coefficients, 
the  standard  deviations  being  lower  than  those  obtained  from  the  eccentricity  fit  alone. 
The  values  were: 


oJ.° 

10  C ,  5  = 

-3+13 

9-',° 

10  S. 5  “ 

-4+11 

0.-1, 2 
.0  c,5  = 

-  63  +  15 

0.-1, 2 
,0S>5  = 

-  46  +  18 

These  values  were  used  in  the  solutions  for  the  individual  coefficients,  but  it  was  found 

-l,° 

necessary  to  double  the  standard  deviations  of  the  two  smaller  coefficients,  C.,  and 
_ I  ,0 

.  The  values  of  eccentricity  are  plotted  in  Fig  14  and  the  curve  shows  the  SIMRES 
fitting  to  the  values. 

9_0> ' 

The  values  of  the  (1,0)  terms  from  the  SIMRES  fitting  were  10  C  .  =  -20.7  t  9.2 

9_0, 1 

and  10  Sjj  =  3.2  ±  4.5  ,  with  e  *=  0.960  as  compared  with  0.675  from  the  solution  for 
i  alone.  These  values  were  not  used  because  of  the  higher  value  of  e  and  also  because 
they  did  not  fit  the  solutions  so  well. 

3.18  Ariel  3,  1967-42A  (i  =  80.17°,  e  =  0.007) 

Ariel  3  was  the  first  satellite  to  be  used  for  evaluating  lumped  15th-order 

g 

harmonics,  when  Gooding  analysed  the  281  orbits  he  had  determined  from  Minitrack 
observations.  Subsequently  the  orbit  has  served  as  a  standard  for  testing  in  Gooding’s 
development  of  the  THROE  and  SIMRES  programs.  The  best  available  values’’  are  from  a 
SIMRES  fitting  of  inclination  and  eccentricity  with  ten  pairs  of  coefficients,  which 
gives  the  following  values  for  the  lumped  harmonic  coefficients: 


i 


20 


9-0, 1 

q  0,1 

10  C, 5 

=  -  23.1  ±  1.6 

,o's15 

=  -  8.6  i  1.3 

oJ.O 

oJ.O 

10  C15 

=  -  54.7  ±  3.2 

1°  S15 

=  -  37.2  +  2.6 

q_->,2 

9_-',2 

1°  C15 

=  -  130.6  ±  10.7 

10  S15 

=  -  96.8  t  9.0  . 

,2 

.-1,2 

and 

coefficients 

required  their 

standard  deviations  increased 

1,0  _-l,2 

The  C)5  ,  C)5 

by  a  factor  of  2  in  the  solutions  for  individual  harmonics. 


It  is  of  interest  to  look  back  at  the  values  of 

8 


9- 

10  C 


0,1  9_0,1 

and  10  S  originally 
1  o  1  o 


obtained  by  Gooding  10  years  ago  with  (y,q)  =  (1,0)  only;  they  were  -19.9  ±  1.2  and 
-7.7  ±  0.8.  These  are  not  far  from  the  values  required  by  our  final  solutions  here, 
namely  -22.9  and  -8.0. 


3.19  Meteor  3,  1970-19A  (i  »  81.16  ,  e  =  0.005) 


Meteor  3  passed  through  exact  15th-order  resonance  on  4  July  1979,  and  the  changes  in 
both  inclination  and  eccentricity  have  been  analysed  using  34  US  Navy  orbits  over  a  period 
from  II  March  to  28  October  1979.  During  this  time  £  changed  from  -7  to  +8  deg/day. 


The  values  of  inclination,  cleared  of  perturbations  except  those  due  to  resonance, 
were  first  fitted  with  THROE  using  just  (y,q)  =  (1,0).  This  seemed  a  fairly  satisfactory 
fit,  with  e  =  0.413.  However,  another  THROE  run  with  (y,q)  =  (1,0),  (1,1)  and  (1,-1) 
was  tried,  and,  although  the  value  of  eccentricity  was  only  0.005,  the  e-terms  in 
equation  (8)  were  determinate,  and  E  was  reduced  to  0.297.  The  (1,0)  terms  from  this 
second  fit,  with  lumped  coefficients 
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0,1 

15 


21.0  ±  1.6 


9 

10  S 


0,1 

15 


-  I  .  I  ±  1 .3  , 


were  accented,  as  they  gave  better  results  in  the  solutions  for  the  individual 

_°,1 

coefficients,  although  the  standard  deviation  of  the  value  still  had  to  be 

increased  by  a  factor  of  4.  The  values  of  inclination  are  plotted  in  Fig  15  and  the 
curve  shows  the  THROE  fitting  to  the  values. 


The  values  of  eccentricity  were  fitted  with  THROE  over  the  same  period  with  (v,q)  = 
(1,1)  and  (1,-1).  The  inclusion  of  (0,1)  terms  gave  no  advantage.  A  further  fitting 
with  (y,q)  ■*  (1,0),  (1,1)  and  (1,-1),  was  used  with  the  corresponding  THROE  fitting  of 
inclination  in  a  SIMRES  fit.  This  yielded  values  of  the  lumped  coefficients  with 
standard  deviations  less  than  those  from  the  THROE  fitting  with  (v,q)  =  (1,1)  and  (1,-1). 
The  values  were: 


o_ 1  >0 

1°  0.5  " 

-  26.2  ±  4.1 

9-1 
10  S 

1,0 

15 

9-->,2 

9-" 

■1,2 

1°  0.5  - 

-  128  ±  29 

10  Sj 

15 

15.0  +  5.0 


130  i  37 
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These  values  were  used  in  the  individual  solutions  but  the  standard  deviations  of  C 
_1,0  1 
and  had  to  be  increased  by  factors  of  10  and  4  respectively.  The  values  of 

eccentricity  and  the  curve  given  by  the  SIMRES  fit  are  shown  in  Fig  16. 

_0,1  _0,1 

The  values  of  the  and  lumped  coefficients  from  the  SIMRES  fitting 

were  not  used  because  they  had  larger  standard  deviations  than  the  values  from  the  THROE 
fitting  of  inclination  alone  and  they  did  not  fit  the  solutions  so  well. 

3.20  OGO  4,  1967-73A  (i  =  85.98°,  e  =  0.025) 

The  Orbiting  Geophysical  Observatory  4  passed  through  15th-order  resonance  on 

2 

22  May  1970  and  the  values  of  inclination  from  27  US  Navy  orbits  were  analysed  in  1975 
The  values  obtained  then  from  a  THROE  fitting  with  (y,q)  =  (1,0)  and  (1,1),  the  (1,-1) 
terms  being  indeterminate,  has  been  accepted,  the  values  being: 

o_0,l  q_0, 1 

10  C)5  =  -  13.9  ±  2.3  10  S15  =  “  6-4  1  3-3  • 


In  the  previous  evaluation  ,  the  standard  deviations  had  to  be  relaxed;  but  here  they 
were  used  unchanged. 

The  values  of  eccentricity  were  not  analysed  for  the  previous  determination  of  even 
harmonic  coefficients',  but  they  have  now  been  utilised.  A  fitting  with  THROE  using 
(y,q)  =  0.1)  and  (1,-1)  was  first  tried  with  little  success,  c  being  3.57.  However, 
when  a  (0,1)  term  was  included,  for  the  reason  given  in  section  3.1,  e  was  reduced  to 
1.47.  The  values  for  the  lumped  coefficients  are  as  follows: 


9- 
10  C 


1,0 

15 


9--1 *2 
10  C>5 


=  -  87  t  38 


122  ±  65 


oJ.O 

10*S]5  =  84+60 


10  S)5  =  -  1 75  ±  86 


These  values  were  used  in  the  solutions  for  the  individual  harmonics  but  it  was  found 

J,0 

necessary  to  increase  the  standard  deviation  of  by  a  factor  of  2.  The  values  of 

eccentricity  and  the  fitting  by  THROE  are  shown  in  Fig  17. 

3.21  SESP  I,  I971-54A  (i  =  90.21°,  e  =  0.002) 

This  orbit  passed  very  slowly  through  15th-order  resonance  between  1972  and  1977, 
and  King-Hele*’  analysed  the  variations  in  inclination  and  eccentricity  from  269  weekly 
US  Navy  orbits  to  give  the  following  values  of  lumped  coefficients: 


nJV 

1°  C15  " 

-  16.40  ±  0.24 

9-°.' 

10  S,5  = 

-  5.37  +  0.15 

9J.° 

10C,5  = 

-92+48 

qJ>0 

.0  S , 5  = 

-  170  +  56 

9--'. 2 

1°  C,5  ’ 

-  62.9  +  2.6 

9--'. 2 

,0  S15  - 

-  53.4  +  1.6 

The  first  two  coefficients  came  from  a  fitting  with  (y,q)  -  (1,0)  and  (2,0). 
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The  two  eccentricity  resonances  were  four  years  apart  and  were  analysed  separately. 

The  resonance  associated  with  (y,q)  =  (1,1)  was  very  weak  because  F,,  „  =  0  at 

o  10  1  o  j  I  j  j  o 

i  =  90  ,  and  consequently  the  (C,S)  ’  coefficients  are  of  very  poor  accuracy:  the 

J.O 

standard  deviation  of  had  to  be  doubled  in  the  solutions  for  individual  coeffi¬ 

cients.  The  last  two  values  came  from  a  fitting  with  (y,q)  =  (1,-1),  (2,-1)  and  (1,0). 

3.22  Nimbus  1  rocket,  I964-52B  (i  =  98.68°,  e  =  0.023) 

The  orbit  of  Nimbus  1  rocket  passed  through  l5th-order  resonance  on  5  June  1970, 
and  more  than  2000  observations  were  used  by  Hiller'2  to  determine  the  orbit  with  PROP 
at  25  epochs  between  March  and  September  1970.  The  25  values  of  inclination,  together 
with  16  from  US  Navy  orbits, were  fitted  using  THROE  with  (y,q)  =  (1,0),  (2,0)  and  (1,1) 


to  give: 


28.3  ±  2.0 


=  1.5  +  2.0 


The  variation  in  inclination,  shown  in  Fig  3  of  Ref  17,  was  surprisingly  small.  The  C 
value  fitted  well  in  the  solutions  for  individual  coefficients,  but  the  S  value  did  not, 
and  its  standard  deviation  had  to  be  increased  by  a  factor  of  4. 

Hiller  also  fitted  the  values  of  eccentricity  and  inclination  together  by  SIMRES 
with  (y,q)  =  (1,0),  (2,0),  (1,1)  and  (1,-1),  and  he  obtained  the  following  values  of 
lumped  harmonics: 


qj.0 

,0C15  = 

-  88  i  7 

oJ>° 

10S,5  * 

-37  +  8 

9-_1  ’2 

1°  C15  * 

-3  +  5 

q_-l,2 

1°  s15  = 

-  34  ±  11 

(These  values  are  obtained  from  the  coefficients  -B,  A,  -D  and  C  respectively  derived  in 
Ref  17,  after  dividing  by  0.5141  for  the  first  pair  and  -0.4713  for  the  second  pair.) 

The  standard  deviations  of  the  C  values  had  to  be  increased  by  factors  of  4  and  2 
respectively  when  used  in  the  solutions.  The  fitting  of  e  is  shown  in  Fig  5  of 
Ref  17. 

3.23  0V1-8,  1 966-63A  (i  =  144.16°,  e  =  0,003) 

This  satellite  passed  through  exact  15th-order  resonance  on  19  December  1976.  There 

were  25  US  Navy  orbits  available  for  analysis  from  9  October  1976  to  27  March  1977  and 

during  this  period  $  changed  from  -II  to  +14  deg/day. 

The  25  values  of  inclination,  cleared  of  perturbations  except  those  due  to 
resonance,  were  fitted  by  THROE  using  (y,q)  =  (1,0).  This  fitting  gave  e  =  0.470  and  the 
values  of  the  lumped  harmonics  were: 

9-°.'  oJM 

10  C|5  =  72000  ‘  16800  10  S  =  -  5960  ±  7630  . 

_°» 1 

The  standard  deviation  on  the  Sjj  value  had  to  be  increased  by  a  factor  of  4  in  the 

individual  solutions.  The  values  of  inclination  and  the  curve  showing  the  THROE  fit  with 

(>,q)  "  (1,0)  are  given  in  Fig  18. 
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An  analysis  of  the  eccentricity  was  not  attempted  as  the  values  were  not  considered 
to  be  accurate  enough. 

3.24  Orbits  not  used 

2 

Two  of  the  orbits  used  in  our  previous  determination  ,  namely  1964-05A  and  1970-65D 
have  been  discarded.  The  first,  at  inclination  31.5°,  is  superseded  by  I965-09A,  a  much 
better  specimen  at  virtually  the  same  inclination.  The  second,  1970-65D  at  51.2°  inclina 
tion,  was  a  high-drag  orbit  included  previously  only  because  of  a  dearth  of  orbits  at 
inclinations  near  50°.  Now  that  1974-34A  is  available,  there  is  no  need  for  1970-65D. 

Two  new  orbits,  1973-99A  and  1975-67A,  both  at  inclination  71.0°,  were  analysed  in 
the  hope  of  improving  the  coverage  at  inclinations  between  65.7°  and  74.0°,  where  we  have 
only  one  (high-drag)  orbit. 

The  analysis  of  1973-99A  utilised  20  US  Navy  values  of  inclination  between 
11  February  and  23  June  1974  (resonance  was  on  29  April);  but  no  clear  pattern  of  varia¬ 
tion  emerged,  and  the  values  of  lumped  coefficients  were  indeterminate. 

The  analysis  of  1975-67A  was  based  on  16  US  Navy  values  between  21  September  1975 
and  4  January  1976  (resonance  being  on  7  October  1975).  The  results  were  slightly  better 
than  for  I973-99A,  but  the  variation  in  inclination  was  very  feeble  and  the  values  of  the 
two  lumped  harmonics  in  a  (y,q)  =  (1,0)  fitting  were: 

9-0,1  9-0’1 

10  C(5  =  -3+10  IOyS)5  =  -5+5  , 

with  e  =  0.62.  These  values  were  not  considered  accurate  enough  to  be  worth  using, 
although  in  retrospect  it  is  apparent  that  the  S  value  would  have  fitted  well. 

1 8 

We  decided  not  to  use  the  values  obtained  by  Wagner  and  Klosko  from  1971-83B 

2 

(included  in  our  previous  determination  ),  partly  because  we  had  already  had  a  good 
satellite  at  a  similar  inclination  (33°)  and  partly  because  we  wished  to  avoid  any 
direct  link  with  the  Goddard  Earth  Models.  In  assessing  the  accuracy  of  Earth  models, 
it  is  essential  to  try  to  make  independent  evaluations. 

For  the  same  reasons  we  decided  not  to  use  the  values  obtained  from  the  analyses  of 
inclination  at  15th-order  resonance  by  Klokotnik 1 ^  These  are  all  at  inclinations 

close  to  those  already  represented  in  our  analyses  and  are  most  useful  as  an  independent 
check  on  our  procedures. 

4  THE  SOLUTIONS  FOR  INDIVIDUAL  COEFFICIENTS 
4 . I  The  equations  to  be  solved 

Each  of  the  lumped  coefficients  derived  in  sections  3.1  to  3.23  can  be  expressed 
as  a  linear  sum  of  individual  l5th-order  coefficients,  C,  and  S.  ,  by  equations 
(9)  to  (II),  which,  on  reverting  to  the  Q-notation  of  equation  (5),  may  more  compactly 
be  written: 
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with  similar  equations  for  the  S  coefficients.  The  Q  constants  have  been  evaluated 
with  the  computer  program  PROF,  and  all  the  relevant  values  for  the  satellites  used  are 
given  in  Tables  3  to  5  (pages  33  to  35).  There  are  23  satellites  yielding  values  of  the 
form  (13)  for  odd-degree  harmonics,  and  16  of  them  also  give  equations  of  the  form  (14), 
leading  to  32  equations  for  even-degree  harmonics. 


Following  the  procedure  which  proved  successful  previously,  we  add  constraint 
equations  of  the  form 


Et,l 5  *  0  4  10"5/l2 

SM5  -  0  *  hTW 


(15) 


where  8.=  15,  17,  19  ...  for  the  odd-degree  harmonics,  and  i  =  16,  18,  20  ...  for 

3 

those  of  even  degree.  These  equations  express  the  expectation  that  the  order  of 

-5  2 

magnitude  of  the  individual  coefficients  of  degree  S.  is  10  )l  for  15  <  i  <  50  , 
as  is  confirmed  in  a  general  way  by  the  Goddard  Earth  Model  10C  (Ref  22). 

Thus,  when  solving  for  N  harmonics,  we  have  23  +  N  pairs  of  equations  for  odd- 
degree  harmonics,  and  32  +  N  pairs  of  equations  for  even-degree  harmonics. 

4.2  The  method  of  solution  -  a  modified  least-squares 


Our  23  satellites  give  results  of  immensely  variable  accuracy  and  reliability. 

Some,  in  particular  1974-34A,  1963-24B,  1970-1UA,  1971-1 3B  and  1971-54A,  are  of  low  drag 
and  give  accurate  results.  Some,  such  as  1964-84A,  1977-12B,  1971-106A  and  1971-18B,  are 
of  high  drag  and  are  included  because  they  are  the  only  satellites  available  to  fill  gaps 
in  the  coverage  of  inclination.  Other  orbits  fall  between  these  extremes. 


When  the  equations  were  first  solved  by  least  squares,  it  was  found,  as  expected, 
that  some  of  the  values  from  the  less  reliable  satellites  did  not  fit.  Of  the  64  values 
of  the  coefficients  with  (q,k)  =  (1,0)  and  (-1,2),  three  appeared  to  be  in  error  by 
about  10  sd  and  another  seven  by  nearly  5  sd . 


To  attempt  a  straightforward  least-squares  solution  is  an  inadequate  response  to 
such  an  abnormal  distribution,  because  large  spurious  values  have  too  much  power  in  a 
least-squares  fit;  or,  to  put  it  more  ecologically,  the  normal  distribution  is  the  proper 
habitat  for  least-squares  procedures.  One  possibility  would  be  to  reject  at  an  arbitrary 
level;  but  this  is  a  drastic  and  unsophisticated  procedure,  a  blunt  instrument  which  not 
only  produces  discontinuities  as  the  rejection  level  varies,  but  also  could  easily 
nullify  our  efforts  to  provide  some  representation  for  unfashionable  inclinations.  So  we 
rejected  all-or-none  rejection,  and  instead  introduced  two  quantum  jumps,  increasing  the 
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standard  deviations  of  ill-fitting  values  by  factors  of  either  2  or  4.  (We  also  allowed 
an  increase  by  a  factor  of  10,  which  is  closely  equivalent  to  complete  rejection  and 
needs  no  further  comment.)  This  process  ensures  that,  while  no  ill-fitting  value  has  too 
much  power,  the  constraint  implied  by  its  presence  is  still  operative,  though  weaker. 

Also  the  distribution  is  brought  much  closer  to  normality,  thus  giving  the  least-squares 
process  a  more  natural  habitat. 

After  we  had  computed  and  considered  a  large  number  of  solutions  of  the  equations, 
it  became  clear  that  the  optimum  number  of  coefficients  would  be  between  10  and  13  for 
both  odd-  and  even-degree  C  and  S  coefficients.  In  these  circumstances  an  acceptably 
distributed  set  of  residuals  (with  one  exception)  was  obtained  by  applying  the  quantum 
multipliers  (2  or  4)  so  as  to  keep  the  weighted  residuals  from  each  lumped  harmonic  less 
than  1.5.  With  this  choice,  the  values  of  the  measure  of  fit  •  were  all  between  0.8 
and  1.0,  and  all  but  one  of  the  individual  C  and  S  coefficients  had  weighted  residual 
less  than  1.5.  Of  the  46  odd-degree  lumped  harmonics,  three  needed  their  sd  multiplied 
by  a  factor  of  2,  and  six  needed  a  factor  of  4.  Of  the  64  even-degree  lumped  coefficient 
twelve  needed  a  factor  of  2,  seven  needed  a  factor  of  4  and  three  had  to  be  multiplied 
by  10. 

4 . 3  The  solutions  for  individual  coefficients  of  odd  degree 

When  the  23  equations  of  type  (13)  and  N  equations  of  type  (15)  were  solved  by- 
least  squares  for  N  coefficients,  the  values  of  the  measure  of  fit  e  for  7  ^  N  <  1  t 
we  re : 


N 

7 

8 

9 

10 

1  1 

12 

13 

c 

equations 

3.91 

2.27 

1.15 

0.99 

0.93 

0.92 

0.92 

s 

equations 

1.21 

1.16 

1  .09 

0.84 

0.83 

0.82 

0.82 

As  before. 

k"  is 

the  sum  of 

squares 

of  weighted 

residuals 

divided 

by  the  number  of 

degrees  of  freedom,  and  the  weighted  residual  is  the  residual  for  each  lumped  coefficient 
divided  by  the  standard  deviation  for  that  coefficient  as  given  in  Table  1. 

It  is  obviously  advantageous,  for  both  C  and  S  ,  to  solve  for  at  least  ten 
coefficients:  it  is  also  necessary,  because  the  Q  factors  in  equation  (13)  remain  quite 
large  up  to  the  10th  for  the  satellites  of  lowest  inclination,  as  Table  3  shows.  Til's  i> 
confirmed  by  the  behaviour  of  the  solution:  in  the  9-coefficient  S  solution  the 
weighted  residuals  for  the  first  four  satellites  in  Table  1  all  numerically  exceed  1.0. 
but  in  the  10-coefficient  solutions  they  are  all  less  than  0.7. 

The  solutions  for  13  coefficients  offer  no  advantage  over  those  for  12  coefficients 
so  the  choice  lies  between  the  10-,  11-  and  12-coefficient  solutions.  The  1 1 -coef f i cient 
solution  seems  preferable  to  the  10-coefficient,  because  of  the  6%  decrease  in  >:  for 
the  C  equations;  but  the  choice  between  the  11-  and  12-coefficient  solutions  is 

difficult.  The  onlv  point  in  favour  of  12  coefficients  is  that,  for  one  satellite, 

°,1- 

1  uqq-n 8 B ,  the  value  of  the  12th  term  on  the  right  hand  side  of  (13),  ^  , 

exceeds  the  standard  deviation  allocated  to  this  satellite,  and  might  therefore  seem  to 
be  needed  -  but  for  the  fact  that  the  I  I -coef f i cient  solution  fits  well.  The  arguments 
against  12  coefficients  are  Lhat  t  decreases  very  little  and  that  the  12th  coefficient 
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is  formally  indeterminate  for  both  C  and  S  .  Since  only  one  value  (that  of  C„,  ) 

Z  /  ,  1  j 

changes  by  more  than  0.2  sd  on  going  from  li  to  12  coefficients,  the  choice  is  not 
crucial,  and,  as  it  is  not  crucial,  the  lower  number  of  coefficients  is  to  be  preferred. 
The  1 1 -coefficient  solutions  are  given  in  Table  6. 


Table  6 


The  values  of  odd-degree  ^  and  S^  ^  given  by  the  1 1 -coefficient  solutions 


e 

l0\l5 

l0\l5 

15 

-22.7  ±  0.6 

-7.4  t  0.6 

17 

11.3  ±  1.0 

6.7  t  1.2 

19 

-13.3  ±  0.8 

-11.8  t  0.9 

21 

15.9  ±  0.7 

8.7  t  0.8 

23 

14.3  ±  1.6 

-1.3  i  1.9 

25 

-12.7  i  2.0 

0.6  +  2.4 

27 

-6.8  t  1.4 

12.7  *  2.0 

29 

00 

+1 

<M 

CM 

1 

0.3  .*  1.9 

31 

27.9  ±  2.9 

-2.0  i  3.9 

33 

6.5  ±  2.9 

-12.0  *  3.8 

35 

-6.8  ±  4.1 

3.3  +  4.6 

The  weighted  residuals  in  the  23  satellite  equations  (13)  and  the  II  constraint 
equations  (15)  are  given  in  Table  7.  All  the  residuals  of  the  lumped  harmonics  are  less 
than  1.5  as  a  result  of  applying  the  quantum  multipliers,  as  explained  in  section  4.2. 


Table  7 

Weighted  residuals  in  the  34  equations  for  odd-degree  harmonics, 
from  the  11-coefficient  solutions 


Satellite  equations 


Constraint 

equations 

Degree  E  of 

coefficient 

S.,15 

se,!5 

15 

0.51 

0.17 

!  17 

-0.33 

-0.19 

19 

0.48 

0.42 

j  21 

-0.70 

-0.38 

23 

-0.75 

0.07 

25 

0.79 

-0.04 

27 

0.49 

-0.93 

|  29 

0.19 

-0.03 

31 

-2.68 

0.19 

33 

-0.70 

1.31 

35 

0.83 

-0.41 

Satellite 

0,1 

C 1 5 

_0,l 

S  1  5 

65-09A 

-0.  15 

0.06 

69-68B 

0.02 

0.06 

64-84A 

-0.17 

-0.78 

79-82A 

0.03 

0.38 

71-30B 

-0.92 

-1  .01 

74-34A 

0.  1  1 

-0.  14 

71-58B 

0.85 

0.85 

62-I5A 

-0.40 

0.18 

65-53B 

-0.05 

-1.27 

63-24B 

-0.04 

0.34 

70-87A 

0.55 

-0.33 

77-J2B 

-1.17 

-0.07 

7 1  -1 06 A 

-0.81 

1  .45 

71-I8B 

-0.80 

0.94 

70-1 1  1A 

-0.12 

0.24 

71-13B 

0.91 

-0.51 

77-95B 

0.31 

0.92 

67-42A 

-0.27 

-0.48 

70-1 9A 

0.48 

1.18 

67-73A 

0.32 

-0.50 

71-54A 

-0.07 

0.03 

64-52B 

-1.00 

1.26 

66-6  3A 

1.38 

-0.82 

One  notable  feature  of  Table  7  is  the  excellent  fitting  of  nearly  all  the  accurate 
satellites.  The  weighted  residuals  on  both  the  C  ant  S  equations  are  less  than  0.3' 
for  1 965-09A,  1969-b8B,  1974-34A,  19M-24B,  1970-1 1 1A  and  197I-54A,  and  there  are  no 
relaxations  of  the  standard  deviation  on  any  of  these. 

The  second  outstanding  feature  is  the  very  large  value  of  ^  ,  which  is  also 

apparent  of  course  in  Table  6.  This  large  value  caused  us  much  concern,  and  we 
arbitrarily  altered  the  standard  deviations  of  many  of  the  lumped  harmonics  by  a  factor 
of  10  to  try  to  identify  the  'culprit'  responsible  for  this  high  value;  but  all  the 
satellites  indicted  had  to  be  acquitted,  and  the  high  value  seems  to  result  from  a 

consensus.  We  also  computed  solutions  with  the  constraints  on  C  _  relaxed  to 

“5  2  -  .  ^  ,  I  5 

2  x  10  /i  ;  inevitably  the  value  of  increased  (from  27.9  to  30.0  »  10  ),  the 

value  of  f  was  reduced  (from  0.93  to  0.70),  and  the  standard  deviations  were  also  much 

reduced.  This  solution  was  not  accepted  because  it  was  not  compatible  with  the  S 

solutions.  Since  making  the  constraint  twice  as  stringent  only  reduces  this  large 

coefficient  from  30.0  to  27.9,  we  are  forced  to  the  conclusion  that  the  value  is  realistic 

and  is  likely  to  be  greater,  not  less,  than  27.9. 

The  values  of  the  lumped  harmonics  from  22  of  the  satellites  are  plotted  against 
inclination  in  Fig  19,  after  multiplication  by  F^  ^  ^  to  keep  the  numerical  values 
to  a  reasonable  level.  The  satellite  omitted  is  1966-63A,  because  its  inclination  is 
144°.  The  standard  deviations  in  Table  1  are  marked  as  bars.  The  curves  in  Fig  19  shea¬ 
the  variations  given  by  the  1 l-roefficient  solutions.  The  fitting  is  quite  sat i sf actor\ 
and  it  is  evident  that  some  of  the  less  accurate  values  were  in  need  of  the  increase  in 
standard  deviation. 

4.4  The  solutions  for  individual  coefficients  of  even  degree 


;quat 

ions  of 

type 

(15) 

the 

measure 

of  fi 

.t  € 

10 

1 1 

12 

13 

0.99 

0.99 

0.98 

0.98 

0.91 

0.91 

0.91 

0.91 

The  4%  decrease  in  t  between  the  9-  and  10-coefficient  C  solutions  is 

substantial;  but  there  is  not  much  to  be  said  in  favour  of  more  than  10  coefficients. 

_-',2 

The  value  of  the  11th  0-coefl‘icient  is  most  significant  with  C  for  1971  —5 4 A ,  but 

-1,2_ 

the  value  of  Q  C1£  c  in  the  I  I -coe f f icient  solution  is  only  one-fifth  of  the 
36  3b  ,  I  5 

standard  deviation.  So  the  Nth  coefficient  should  not  be  required.  This  is  confirmed 
by  the  fact  that  the  values  of  the  individual  coefficients  are  not  appreciably  altered 
by  increasing  N  beyond  10.  So  we  choose  the  10-coefficient  solutions,  given  in 
Table  8. 
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Table  8 


The  values  of  even-degree  C  , 

_ 1  J 


and  S 


SL,  15 


given  by  the  10-coefficient  solutions 


i 

,0\,5 

'°95£,I5 

16 

-11.0  ±  2.7 

-21.5  ±  1.7 

18 

-43.0  ±  1.8 

-22.5  ±  1.2 

20 

-24.3  ±  2.3 

-6.2  ±  1.6 

22 

24.1  ±  2.0 

10.2  ±  1.6 

24 

1.4  ±  3.8 

-21.8  ±  3.3 

26 

-13.3  ±  5.8 

14.4  ±  5.5 

28 

-15.4  ±  6.4 

-8.4  ±  6.3 

30 

-4.0  ±  6.8 

-16.0  ±  6.2 

32 

7.8  ±  6.2 

2.5  ±  5.1 

34 

9.6  ±  6.3 

5.6  ±  5.2 

The  weighted  residuals  in  the  32  satellite  equations  (14)  and  the  10  constraint 
equations  (15)  are  given  in  Table  9;  all  the  84  weighted  residuals  are  less  than  1.5. 


Table  9 

Weighted  residuals  in  the  42  equations  for  even-degree  harmonics, 
from  the  10-coefficient  solutions 


Satellite  equations 

Constraint 

equations 

Satellite 

J>° 

C,5 

.-1.2 

C 1 5 

J,o 

S1 5 

_—  1  » 2 

S1 5 

Degree  i  of 
coefficient 

^£,15 

S£,  1  5 

79-82A 

-0.08 

-0.13 

-0.30 

0.85 

16 

0.28 

0.55 

71-30B 

0.80 

1.07 

-1.33 

0.59 

18 

1.40 

0.73 

74-34A 

-0.26 

-0.68 

-0.23 

0.01 

20 

0.97 

0.25 

71-58B 

-1.19 

-l  .38 

0.94 

-0.64 

22 

-1 . 16 

-0.49 

62-I5A 

0.80 

1.24 

0.91 

-0.12 

24 

-0.08 

1.25 

65-53B 

1.33 

0.45 

0.26 

-0.88 

26 

0.90 

-0.97 

63-24B 

-0.97 

0.66 

0.24 

0.53 

28 

1.20 

0.65 

7 1-1 06 A 

-0.80 

-1.24 

-0.95 

0.12 

30 

0.36 

1 .44 

70-1 1 1A 

0.25 

-0.23 

-0.77 

-1.27 

32 

-0.80 

-0.26 

71-13B 

-0.46 

0.34 

0.01 

0.36 

34 

-1.10 

0.65 

77-95B 

0.82 

-0.26 

1 .01 

-0.20 

67-42A 

-1.21 

-1.38 

-0.40 

-0.96 

70-19A 

0.67 

-0.29 

1.27 

-0.82 

67-73A 

-0.28 

-0.85 

1.23 

-1.4! 

7 1 -54A 

-0.23 

0.60 

-0.87 

0.20; 

64-52B 

-1.22 

1.38 

0.50 

-1.14 

J,0  _  _ —  I » 2 

The  values  of  F.,  . .  „C.,  and  F,,  ,C  are  plotted  against  inclination  in 

1  D  i  I  j  )  o  !_)  I  u  }  I  J  j  /  13 

Fig  20,  and  Fig  21  is  a  similar  diagram  for  the  S  coefficients.  It  should  be 

remembered  that  the  two  sets  of  values  in  Fig  20  are  being  fitted  simultaneously:  thus 

1 ,0 

—  ’  #  o 

the  rather  perverse -looking  course  of  the  curve  for  C.,  in  Fig  20  near  the  i  *  65.7 

,2 

point  is  caused  by  the  need  to  fit  one  of  the  C,^  values  (i  =  58.2  )  which  has 
similar  coefficients.  If  we  make  allowance  for  this  effect,  the  fitting  in  Figs  20  and 
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_-l,2 

21  is  entirely  satisfactory,  being  particularly  good.  It  is  also  apparent  that 

a  few  of  the  values  thoroughly  deserve  their  increased  standard  deviation. 

5  DISCUSSION 

In  comparing  our  new  values  for  the  15th-order  coefficients  with  previous  values, 

there  is  a  problem  in  deciding  which  previous  sets  of  coefficients  should  be  considered. 

^23 

Klokocnik  and  Pospisilova  have  pointed  out  that  the  variation  of  lumped  coefficients 

with  inclination  provides  a  useful  test  of  the  accuracy  of  existing  geopotential  models, 

because  the  values  predicted  by  the  models  can  be  compared  with  each  other  and  also  with 

the  more  accurate  values  derived  independently  by  resonance  analysis  for  particular 

inclinations.  Klokocnfk  and  Pospisilova  present  the  variations  of  selected  lumped 

coefficients,  including  ^  ^(C,S)^’'  ,  as  given  by  11  different  geopotential  models 

derived  in  the  past  10  years.  Of  these  11  models,  four  can  be  excluded  because  they  do 

not  go  beyond  degree  18.  Of  the  remaining  seven,  three  are  Goddard  Earth  Models,  of  which 
24 

the  latest,  GEM  10B  ,  can  be  taken  as  superseding  the  earlier  ones.  (We  do  not  require 

22 

GEM  10C  ,  which  is  the  same  as  GEM  I  OB  to  degree  36,  because  we  do  not  go  beyond 

23 

degree  36  in  our  solutions.)  We  also  exclude  HARMOGRAV,  which  appears  to  be  the  least 
accurate  of  the  remaining  models. 

This  process  of  elimination  leaves  us  with  four  models,  namely  GEM  10B;  the 

25  „ 

Smithsonian  Standard  Earth  IV. 3  (Gaposchkin  ,  labelled  as  SE  '5'  by  Klokocnfk  and 

v  26 

Pospisilova);  GRIM  2  (the  European  model  );  and  the  purely  terrestrial  model  'Rape  I  1  " 

_0,1  _  _0,1 

(Ref  27).  The  variations  of  ^  and  F | j  ^  7^15  for  fhese  four  models, 

as  given  by  KlokoJnfk  and  Pospi'SilovS,  are  reproduced  in  Figs  22  and  23.  Our  values  of 

9 —  9  — 

10  FC  and  10  FS  from  Table  1  are  also  plotted  whenever  their  standard  deviation  is 
less  than  1 .0,  and  are  marked  as  black  circles  of  diameter  1.0.  Fig  23  shows  that  the 
values  of  FS  given  by  SSE  IV. 3  and  Rapp  1977  cannot  be  regarded  as  realistic,  because 
they  differ  so  greatly  from  those  experienced  by  the  resonant  orbits  for  inclinations 
between  70°  and  90°.  At  90°,  for  example,  1971 —5 A A  gives  lO^FS  =  -3.1  t  0.1  ,  while 
SSE  IV. 3  gives  +24  and  Rapp  1977  gives  +18. 

So  we  are  reduced  to  GEM  I 0B  and  GRIM  2  for  our  comparisons..  Since  GRIM  2 
utilized  our  previous  values  for  15th-order  harmonics,  almost  unchanged,  it  does  not 
provide  an  independent  test.  But  it  does  show  approximately  how  a  curve  based  on  our 
old  values  of  ^  and  S?  ^  fits  the  new  data.  Looking  first  at  Fig  23  (which  is 

the  worse),  we  see  that  the  first  two  points  conform  with  the  dot-dash  curve,  but  the 
three  new  points  at  inclinations  of  43.6°,  50.6°  and  58.2°  do  not;  there  are  smaller 
differences  for  some  of  the  80°-90°  points  where  we  have  used  new  values.  In  Fig  22  the 
situation  is  similar,  except  that  the  58°  point  fits  well.  All  this  is  as  expected  and 
merely  confirms  the  truism  that  new  data  demand  a  new  solution. 

The  only  useful  comparison  therefore  is  with  GEM  10B,  which  is  believed  to  be 
independent  of  our  results.  Figs  22  and  23  show  that  GEM  I0B  agrees  quite  well  with  the 
values  from  the  resonant  orbits  and  can  therefore  be  accepted  as  a  fairly  good  representa¬ 
tion  of  the  llth-order  harmonics  of  odd  degree.  The  individual  odd  degree  values  in 
GEM  I0B  differ  considerably  from  ours,  however,  for  degree  25  and  higher,  and  only  time 
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will  tell  which  is  the  better  set  of  individual  coefficients;  the  average  standard 

-9 

deviation  of  our  odd-degree  coefficients  for  degrees  25  to  35  is  2.8  *  10  ,  while  the 

standard  deviation  of  the  GEM  10B  values  is  estimated^’^  as  3  *  10  9.  For  the  odd-degree 
harmonics  of  degree  15  to  23,  the  values  of  the  coefficients  can  be  regarded  as  quite 
well  established,  and  Table  10  compares  the  values  from  GEM  10B  with  our  solution. 


Table  10 

Comparison  of  odd-degree  15th-order  harmonics  up  to  degree  23, 
given  by  GEM  10B  and  Table  6 


l 

10 

9C 

1,15 

10 

9S 

b2,l5 

GEM  10B 

Table  6 

GEM  10B 

Table  6 

15 

-19.7 

-22.7  ±  0.6 

-6.4 

-7.4  ±  0.6 

17 

2.5 

11.3  ±  1.0 

4.8 

6.7  ±  1.2 

19 

-20.6 

-13.3  ±  0.8 

-15.3 

-11.8  ±  0.9 

21 

15.2 

15.9  ±  0.7 

9.5 

8.7  ±  0.8 

23 

15.4 

14.3  ±  1.6 

4.1 

-1.3  ±  1.9 

.  -9 

If  the  GEM  10B  standard  deviation  is  taken  as  3  *  10  for  all  values,  the  difference 
between  GEM  10B  and  the  corresponding  value  in  our  solution  is,  on  average,  0.8  *  (the 
sum  of  the  two  standard  deviations),  the  S  coefficients  being  in  better  agreement  than 
the  C  coefficients.  If  the  sets  of  values  are  independent,  as  we  believe,  the  agree¬ 
ment  is  very  satisfactory. 

The  variations  of  the  even-degree  GEM  10B  lumped  harmonics  (C,S)]’^  and 

---12..  -  15 
(C,S)jg’  ,  multiplied  by  the  appropriate  F  factors,  are  shown  in  Fig  24.  Also  plotted 

are  those  values  of  lumped  harmonics  determined  from  the  resonant  orbits  (Table  2)  which 

have  standard  deviations  less  chan  2.5.  The  points  are  plotted  as  circles  o£  diameter 

2.5.  Again  the  impression  is  that  GEM  10B  provides  quite  a  good  approximation  to  the 

-1,°  Q 

observed  values,  the  worst  discrepancy  being  for  at  44  inclination. 

The  comparison  of  the  first  five  individual  coefficients,  Table  II,  shows  excellent 
agreement  for  the  C  but  greater  discrepancies  for  the  S  coefficients.  If  the 
standard  deviation  of  the  GEM  I0B  values  is  taken  as  3.0,  the  difference  between  GEM  10B 
and  the  corresponding  value  in  our  solution  is,  on  average,  1.2  x  (the  sum  of  the  two 
standard  deviations)  which  is  rather  high  but  acceptable.  In  view  of  the  good  fitting  of 
the  GEM  10B  curve  to  out  lumped  S  values  in  Fig  24,  it  may  seem  surprising  that  there 
are  such  large  differences  between  the  GEM  10B  ^  values  and  ours.  But  this  is  to 

be  expected  because  our  fitted  curves  (Fig  21)  differ  considerably  from  the  GEM  10B 
curve 8. 
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Table  I  1 

Comparison  of  even-degree  l5th-order  harmonics  up  to  degree  24, 
given  by  GEM  10B  and  Table  8 


l 

10 

\.5 

10 

\.5 

GEM  10B 

Table  8 

GEM  10B 

Table  8 

16 

-14.4 

-11.0  ±  2.7 

-27.8 

-21.5  +  1.7 

18 

-48.3 

-43.0  ±  1.8 

-18.6 

-22.5  t  1.2 

20 

-23.9 

-24.3  ±  2.3 

4.8 

-6.2  t  1.6 

22 

24.1 

24. 1  ±  2.0 

-1.3 

10.2  +  1.6 

24 

3.1 

1.4  ±  3.8 

-5.1 

-21 .8  i  3.3 

Our  new  values  of  the  individual  coefficients  are  formally  much  more  accurate  than 
1  2 

our  previous  values  ’  (the  standard  deviation  being  halved,  on  average);  the  new  values 
should  also  be  more  reliable  because  the  coverage  of  inclination  is  much  better. 

6  CONCLUSIONS 

The  aim  of  this  paper  has  been  to  derive  the  best  possible  values  of  15th-order 
harmonic  coefficients  of  odd  and  even  degree,  solely  from  analysis  of  orbits  which  have 
passed  through  15th-order  resonance.  We  have  tried  to  obtain  as  many  reliable  values  as 
possible  of  lumped  15th-order  harmonics  from  orbits  over  a  wide  range  of  inclinations. 
Useful  values  of  lumped  harmonics  of  odd  degree  were  determined  from  analysis  of  the 
changes  in  inclination  of  22  orbits  at  inclinations  between  30°  and  100°,  and  one  at  144°. 
The  distribution  in  inclination  was  far  more  satisfactory  than  has  been  achieved  before, 
but  there  are  still  some  gaps  that  need  filling.  Useful  values  of  lumped  harmonics  of 
even  degree  were  obtained  by  analysing  the  changes  in  orbital  eccentricity  of  16  of  the 
23  orbits,  and  this  yielded  32  pairs  of  equations  for  coefficients  of  even  degree. 

These  values  of  lumped  harmonics  have  been  used  to  solve  for  individual  coefficients, 
and  the  solutions  chosen  as  the  best,  given  in  Tables  6  and  8,  are  the  1  I -coef f icient 
solution  for  odd  degree  and  the  10-coefficient  for  even  degree,  so  that  the  solution  is 
complete  for  degree  15-35.  The  coefficients  of  degree  15-23  should  be  more  accurate  than 
any  obtained  previously,  but  the  errors  inevitably  increase  as  the  degree  increases 
beyond  24. 

The  values  of  lumped  harmonics  also  provide  an  independent  test  of  comprehensive 
Earth  models.  Figs  22  and  23  suggest  that  only  the  Goddard  Earth  Model  10B  is  useful 
for  comparison.  Tables  10  and  II  compare  our  values  of  15th-order  coefficients  with  those 

of  GEM  I0B  for  degree  up  to  24.  [f  it  is  assumed  that  the  GEM  10B  coefficients  have 

-9 

standard  deviations  of  3  *  10  ,  the  difference  between  the  twenty  GEM  I0B  values  and  the 

corresponding  values  in  our  solution  is,  on  average,  1.0  x  (the  sum  of  the  two  standard 
deviations).  If  the  sets  of  values  are  independent,  as  we  believe,  this  agreement  is 
very  satisfactory. 

The  average  standard  deviation  of  our  coefficients  of  degree  15,  16,  17  ....  23  is 
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VALUES  OF  Q',  QT!  ...  Q,’  FOR  THE  23  SATELLITES 
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Fig  1  1965-09A:  variation  of  inclination  near  15th-order  resonance 
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Fig  3  1964-84A:  variation  of  inclination  near  15th-order  resonance 


Fig  4  1979-82A:  variation  of  inclination  near  15th-order  resonance 
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Fig  5  1979-82A:  variation  of  eccentricity  near  15th-order  resonance 


from  US  Navy  orbits 


Fig  6  1971-30B:  variation  of  inclination  near  15th-order  resonance 
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Fig  7  1971-30B:  variation  of  eccentricity  near  15th-order  resonance 


from  US  Navy  orbits 


Fig  8  1871-58B:  variation  of  inclination  near  15th-order  resonance 
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Fig  9  1971-58B:  variation  of  eccentricity  near  15th-order  resonance 
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Fig  10  1965-53B:  variation  of  eccentricity  near  15th-order  resonance 
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Fig  11  1970-87A:  variation  of  inclination  after  15th-order  resonance 
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Fig  13  1977-95B:  variation  of  inclination  near  15th-order 


Fig  14  1977-95B:  variation  of  eccentricity  near  15th-order  resonance 
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Fig  16  1970-19A:  variation  of  eccentricity  near  15th-order  resonance 


Fig  17  1967-73A:  variation  of  eccentricity  near  15th-order  resonance 


1966 -63 A:  variation  of  inclination  near  15th-order  resonance 
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Fig  24  Variation  with  inclination  of  the  four  even-degree  lumped  harmonics  as  given  by 
GEM  10B,  and  values  from  Table  2 
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